Magnetism of one-dimensional Wigner lattices and its impact on charge order 
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The magnetic phase diagram of the quarter-filled generalized Wigner lattice with nearest- and 
next-nearest-neighbor hopping t\ and t2 is explored. We find a region at negative t2 with fully 
saturated ferromagnetic ground states that we attribute to kinetic exchange. Such interaction 
disfavors antiferromagnetism at f2 < and stems from virtual excitations across the charge gap of 
the Wigner lattice, which is much smaller than the Mott-Hubbard gap oc U . Remarkably, we find a 
strong dependence of the charge structure factor on magnetism even in the limit U ^ oo, in contrast 
to the expectation that charge ordering in the Wigner lattice regime should be well described by 
spinless fermions. Our results, obtained using the density-matrix renormahzation group and exact 
diagonalization, can be transparently explained by means of an effective low-energy Hamiltonian. 
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I. INTRODUCTION 

In a Wigner lattice (WL), long-range Coulomb repul- 
sion dominates over the kinetic energy of electrons and 
leads to strong and well-defined charge order ji Origi- 
nally introduced for the electron gas with a homoge- 
neous neutralizing background, this concept was gener- 
alized to electrons on a lattice by Hubbardi^ Evidence 
for quasi one-dimensional (ID) Wigner lattices has been 
found in organic'^''*'^'® and anorganio ^i^'^'^" chain com- 
pounds, nanowiresii and in carbon nanotubes.— Low- 
dimensional WLs are favored by reduced screening i^ii^ 
Moreover, strong correlations induced by large local Hub- 
bard interaction U suppress screening further and protect 
the long-range nature of the Coulomb repulsion, It is, 
however, not straightforward to distinguish a true WL 
from a quantum-mechanical charge-density wave (CDW) 
simply on the basis of the periodicity of the charge mod- 
ulation. In fact, it turns out that the modulation period 
of the WL coincides with that of the Akp CDW-^iiiiiii^ii^ 
This may be surprising, as the microscopic origin of the 
charge order in the two cases is fundamentally different!^ 
(i) The mechanism for the WL is based solely on the clas- 
sical Coulomb repulsion and is dependent only on the 
charge of electrons and not on their fermionic nature. 
The periodicity follows simply from the configuration of 
charges with minimal energy, (ii) Instead, the quantum 
mechanical CDW depends on the Fermi surface topology 
and the instability and modulation reflects the Fermi mo- 
mentum kp. 

If we allow for nearest-neighbor (NN) and next- 
nearest-neighbor (NNN) hoppings ti and t2, we arrive 
at an even more interesting model which, depending on 
the relative size and sign of ii and ^2, may have an elec- 
tron dispersion with two minima instead of one.^^'"'^^ Ac- 
tually, it has been proposed that the edge-sharing CuO 
chain compounds are described by such models with 



\t2\ > \ti\i^ It is then immediately obvious that in the 
case of a four-Fermi-point topology the periodicities of 
WL and CDW no longer coincide. We note that the 
experimentally observed charge modulations in the edge- 
sharing compounds Nai+a:Cu02— '^'^'^ are strong and 
their periodicity consistent only with that of the WL. 
Another possible type of instability in the presence of 
strong correlations, namely, the 2kp Peierls and spin- 
Peierls modulations,—'^ which requires a distortion of 
the lattice, appears to be ruled out in these systems. 
Doped edge-sharing chains are also building blocks of the 
Ca2+.Y2-,Cu50io^''3 and the Sri4_.Ca,Cu2404i-^'^'24 
systems, and pronounced charge order has been observed 
in these compounds as well. 

Here we shall investigate the intrinsic mechanisms for 
the magnetism of generalized Wigner lattices. It is well 
known that NNN hopping t2 has nontrivial consequences 
for magnetism in the ID Hubbard model at general fill- 
ing and may lead to ferromagnetic (FM) states in cer- 
tain cases j^^iS^'^'^^i^ It should be kept in mind that, 
according to the Lieb-Mattis theorem;^ ferromagnetism 
is excluded at any filling in the ID Hubbard model with 
NN hopping (i.e., t2 = 0). For the ID Hubbard model 
with NN and NNN hopping. Fieri et al^ and Daul 
and Noack^^ found that ferromagnetic ground states ap- 
peared above a critical U in those regions of the ti-t2 
plane where four Fermi points exist. These results were 
obtained in the metallic regime where the relevance of 
Fermi surface topology is suggestive; however, the impli- 
cations for the localized electrons of a Wigner crystal are 
unclear. 

The magnetism of generalized WLs is typically dis- 
cussed in terms of effective Heisenberg models where 
the position of the spins is dictated by the charge or- 
der pattern of the WL.— i^i^ The prevailing superex- 
change interactions are antiferromagnetic. However, 
there are also ferromagnetic couplings in edge-sharing 
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chains due to the Hund interaction at the oxygen hg- 
ands that may be larger than the AF interactions 
and render, e.g., the nearest-neighbor interaction Ji 
ferromagnetic i^^i^i^i These features lead to frustration, 
and the resulting helical spin states have been observed 
in the spin-1/2 edge-sharing chain compounds LiCu202^ 
and NaCuaOa^'^'^s 

In this paper, we show that, for a Wigner crystal at 
quarter-filling, there is another intrinsic mechanism that 
may lead to ferromagnetism. By means of a density- 
matrix renormalization group (DMRG) study of the 
ti—t2 Hubbard- Wigner model, which includes local Hub- 
bard U and long-range Coulomb interactions Vi = V/l, 
we show that there is a regime of fully polarized FM 
states at negative ^2- Subsequently, we derive an effec- 
tive magnetic Hamiltonian for the Wigner-lattice regime, 
i.e. |ti|,|i2| V < U , and show that the emergence 
of ferromagnetism can be explained by an effective ki- 
netic exchange mechanism mediated by NNN hopping t2 ■ 
The associated magnetic exchange constant oc 1^12 /e^ de- 
pends on the sign of t2 and therefore kinetic exchange is 
found to favor ferromagnetism for negative t2 and anti- 
ferromagnetism for positive t2 ■ Kinetic exchange involves 
excitations across the charge gap cq of the WL but not 
across the usually much larger Mott-Hubbard gap ^ U, 
as is the case for AF superexchange or for many realiza- 
tions of FM three-particle ring exchangei^ The charge 
gap eo of the generalized WL depends sensitively on the 
commensurability with the underlying crystalline lattice. 
At quarter-filling this gap is particularly large eo ~ V/2. 

If the charge gap eo ~ of the WL is much larger 
than the hopping amplitudes ti and t2 and any expected 
magnetic couplings, a separation of charge and magnetic 
energy scales appears straightforward. Hence, charge or- 
dering in a WL is usually discussed in terms of spinless 
fermions. Magnetism, e.g., antiferromagnetic (AF) su- 
perexchange or ferromagnetic (FM) Hund's rule^*' and 
three-site ring exchange,"^" is then treated as a pertur- 
bation given a particular charge-ordering pattern. One 
would, however, not expect the magnetic order to have 
a strong impact on the underlying charge order, because 
the magnetic energy scale is so much smaller than the 
dominant Coulomb repulsion for all these processes4i 
The motivation for this work was the initial observation 
that the charge structure of the WL, measured by the 
charge structure factor N{q) at g = tt, is strongly affected 
by electron spin, in disagreement with the calculation for 
spinless fermions. Yet there is a region in the — 12 phase 
diagram at negative t2 where iV(7r) is the same for spin- 
less fermions and for fermions with spin, and, moreover, 
N{tt) does not depend on t2 in that parameter range. 
The obvious conjecture is that the ground state should 
be fully spin-polarized in this regime. 

We show here that, due to the kinetic exchange mech- 
anism, these FM ground states emerge and that the ki- 
netic exchange processes have a surprisingly strong im- 
pact on the charge ordering in spite of the classical origin 
of the WL. Indeed, the AF state at ^2 > has dramati- 
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FIG. 1: (Color online) Schematic depiction of relevant 
oc t\t2/tx) kinetic exchange processes. We start from per- 
fect charge order (a and e), where circles denote empty sites 
(Cu^"'' atoms in Nai+a:Cu02), and arrows denote occupied 
sites (Cu^+). NN hopping ti induces excitations (a— >b and 
e— >f) with two domain walls (dashed lines) and cost eo. Two 
different t2 processes, one with (f^g) and one without (b— »c) 
electron exchange, become possible. For FM spins (triplet 
channel) or spinless fermions, however, these two processes 
(a-^d) and (e^h) cancel exactly because of the relative Fermi 
sign in the next-nearest-neighbor hopping process. 



cally weaker charge order than the ferromagnetic or spin- 
less states. While the charge order is reduced oc t2 for 
positive t2, it does not depend on t2 in the FM regime 
t2 > ^2 > t\ with negative t2 ■ Remarkably, for neg- 
ative t2 values below AF reappears, yet the charge 
order then increases with increasing modulus \t2\. The 
boundaries of the FM phase follow from the effective spin 
Hamiltonian as t"^ ^ —?>t\/U and t\ ^ —{U/eDtf, and 
match the magnetic phase boundaries found using the 
DMRG. This peculiar behavior is due to a purely quan- 
tum effect involving destructive interference of kinetic ex- 
change processes in the FM state due to the Pauli prin- 
ciple and constructive interference for the AF case. 

The paper is organized as follows: After introduc- 
ing the Hubbard- Wigner Hamiltonian in Sec. |TT1 we 
present results for the charge structure factor for spinless 
fermions interacting via long-range Coulomb interaction 
in Sec. IIIII We shall see that the results for spinless 
fermions coincide with results for electrons with spin in 
some region of the ti-t2 phase diagram, while they are 
substantially different in other parts of the phase dia- 
gram. We derive an effective Heisenberg Hamiltonian in 
Sec. llVl and show that the magnetic phase diagram, i.e., 
the appearance of the fully saturated FM phase and its 
phase boundaries, can be naturally explained. Next we 
show that the peculiarities found numerically for the the 
charge structure factor find a straightforward analytical 
description in the framework of the effective Hamiltonian. 
Finally, we discuss and summarize our results in Sec. [V] 
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II. HUBBARD- WIGNER MODEL 

The Hubbard- Wigner Hamiltonian investigated in this 
paper is motivated by the one-dimensional edge-sharing 
CuO-chains>i2ii^ Edge-sharing chains are formed by 
Cu04 squares just as in the CuO planes of the high- 
Tc compounds, but these units are differently linked. 
The edge-sharing arrangement leads to small nearest- 
neighbor hopping matrix element ti due to the almost 
90° Cu-O-Cu coordination, and some contribution to ti 
stems from direct Cu d-d overlap—^ Moreover, the struc- 
ture leads to a comparatively large matrix element t2 
between second neighbor Cu ions stemming from a Cu- 
0-0-Cu path.'^^ Thus edge-sharing chains, in contrast to 
the ISO*' bonded high-T^ cuprates, fulfill the fundamen- 
tal criterion for a WL, namely that the kinetic energy is 
small compared to the nearest-neighbor Coulomb inter- 
action, in an optimal way. While the Coulomb repulsion 
is screened by a static dielectric constant in these insu- 
lators, the one over distance decay of the interaction is 
preserved and must be taken into account. Truncation 
of the interaction may have serious consequences for the 
charge-order pattern of generalized WLs^ as well as for 
their charge excitations 

The relevant states of Cu that need to be included in 
a low-energy model are: Cu'^"'" or, more precisely, the 
Cu d^-ligand hole singlet state, Cu^"*" with spin-1/2, and 
Cu^+, corresponding to the filled d shell. These states 
can be expressed in the frame of a single-orbital Hubbard 
model with 0,1 or 2 electrons per site. Thus we consider 
the Hubbard- Wigner Hamiltonian, where the extension 
Wigner indicates that the long-range Coulomb interac- 
tion is included. This model has the for m^°i"^^ 

H ^ -ti ^(4 ,c,+ i^, + h.c.) - ^2 ^(4,c,+2,. + h.c.) 

L/2 

+ rii^^n^^l + ^ V; ^{rii - n){7ii^i - n) , 

i l—l i 

(1) 

where the operators cj^ (c^ ,j) create (destroy) electrons 
with spin a at lattice site i with i = 1 . . . L. The lo- 
cal density is given by rii^a = c|,o-Cj,o-, "-i = '^ij + "^4' 
and the average density is fi — N^/L for electrons. 
The kinetic energy term includes NN hopping ti and 
NNN hopping ^2, which are both typically much smaller 
than either the on-site Coulomb repulsion U or the NN- 
Coulomb interaction V that parametrizes the long-range 
Coulomb interaction Vi = V/l. In the case of finite rings 
of length L, we truncate the l/l behavior at L/2, which 
is equivalent to replacing V/l by max{V/l, V/{L — l)) for 
< I < L. We have verified that small modifications 
to the l/l behavior do not affect our results. In fact, 
truly long-range Coulomb repulsion is not crucial for the 
results presented here: At quarter filling, both the FM 
kinetic exchange and the weakened charge order can also 
be seen for on-site U and NN Coulomb repulsion Vi only. 



In the following, we consider the model with long-range 
Coulomb interaction and use the NN Coulomb repulsion 
Vi = V as unit of energy. Without loss of generality, ti 
is chosen to be positive. 

The Hamiltonian in Eq.(IT]) contains two ingredients 
that have been shown to favor FM correlations in 
Hubbard-like models: Strong on-site and longer-range 
Coulomb repulsioniS- and, perhaps more importantly, 
NNN hoppingi2^i21i^d& In this paper we address the 
most transparent instance of the WL, namely quarter fill- 
ing n ~ 0.5 for Hamiltonian ([T|). Wc explore the magnetic 
properties within the WL regime, i.e., ti,\t2\ V <^ U, 
which have not been explored before, and find FM ground 
states in a region of the ti-t2 phase diagram with negative 
t2 ■ Quite unexpectedly, we also find a strong influence of 
magnetism on WL charge order. For comparison, we will 
first discuss the charge ordering for spinless fermions at 
half filling, corresponding to the fully spin-polarized case 
with n — 0.5. At small ti and t2, the alternating charge 
order is very rigid and its lowest charge excitations are 
domain walls (DWs) with fractional chargci^'^'^ DWs 
can be induced in a perfectly ordered state via NN hop- 
ping ti, as schematically illustrated in Fig.[TJ Their cre- 
ation costs energy eq ~ V/2 in the case of long-range 
Coulomb interaction, and once created they can move 
easily through the lattice via ti hopping processes. Their 
fractional charge ±1/2 is responsible for the distinctive 
WL features in the optical conductivity and in photoe- 

42 44 4Q 

mission] 1 1^ 

We investigate this model with exact diagonalization 
(ED) for spinless Fermions and chains of up to i = 28 
sites. We use the Lanczos algorithm with a numerical ac- 
curacy of « 10~^ and check its validity by use of full diag- 
onalization for up to 18 sites. These ED calculations were 
done for both the ground state (T = 0) and for small but 
finite temperatures (T = lO^'' — 10"'^), leading to prac- 
tically identical results at small ti and t2 because there 
is a finite charge gap eo ~ V/2, which blocks changes 
at temperatures with an energy scale smaller than this. 
For electrons with spin, we use the DMRG with chains 
of L = 24, 32, 40, and find results consistent with the ED 
for spinless fermions. In the DMRG, we keep 200 to 1200 
states at each step, perform up to 10 finite-size sweeps 
and the neglected weight is < 10~^. For parameters with 
a FM ground state, the energy obtained with the DMRG 
agrees with the ED result. 

III. CHARGE ORDER 

Increasing the NN hopping ti gradually reduces the 
charge ordering^^ until, at ti 0.2V, the charge gap 
vanishes This is refiected in the charge structure factor 

^(9) = {P-qPq)> with pg = l/Ne^exp{-iqr)nr , (2) 

r 

which, for perfect charge alternation, is peaked at q = tt 
with N{tt) = 1. As can be seen in Fig. [21 the results 



4 



1 

0.8 
,0.6 
0.4 
0.2 

















---9-. ? 



0.05 0.1 OAS 0.2 



0.25 



\7 
V 


t2 = -0.1V, U= 100V 




t =-0.1V, U= 4V 

2 


O 


t - +n 1 V II- mnv 

— T'V/. I V f \J — 1 \J\J V 


-1- 


=+0.1 V, U= 4V 




1^=0+0.1 V, spinless 




Eq. (3), FM 




Eq. (3), AF, | = -0.1V 




Eq. (3), AF, | = 0.1V 



FIG. 2: (Color online) Charge structure factor N{q) for q = n 
as a function of nearest-neighbor hopping ti . The dotted line 
for spinless fermions is from ED calculations with L = 28; 
symbols were calculated using the DMRG with L = 24 for 
electrons with spin. Analytic results are obtained from Eq. 
(HI) with A = eo (FM) and A = eo - 2*2 (AF), respectively. 
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FIG. 3: (Color online) The charge structure factor A'^(7r) ver- 
sus t2- Note that it does depend on the sign of t2 for fermions 
with spin and does not for spinless fermions. The results for 
spinless fermions were calculated using exact diagonalization 
{L = 18), and the results for electrons with spin using the 
DMRG {ti = 0.02V, U = 4V, L = 24 and ti = 0.07V, 
U = 100 V, L = 32). 



for spinless fermions, obtained using Lanczos diagonal- 
ization, show that N{7t) is strongly reduced even before 
the gap vanishes, giving a weaker charge density wave. 
We find that, for spinless fermions, the melting of WL 
charge order with ti does not depend on t2 ■ The behav- 
ior of 7V(7r) can be described analytically because only 
few DWs are present at small ti. To leading order, vir- 
tual DW excitations contribute 

^ ~N,2tl/eo (3) 

to the ground state energy. With = L/2, we obtain 



iV(7r) 



1 



1 + (4<i/A)2 ' 



(4) 



given the charge gap A — cq ^ V/2. This expression is 
indicated by the dashed line in Fig. [2] and agrees with 
the numerical data. 

In marked contrast to the gradual change that oc- 
curs with ti, NNN hopping t2 is frustrated for spin- 
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FIG. 4: Phase diagram for spinless fermions determined from 
the charge structure factor N{q), see Eq. ([2]). WL (dark gray): 
strongly charge-ordered WL with A''(7r) > 0.7. tt-CDW (light 
gray): CDW with periodicity tt, but A'^(7r) < 0.7. (This choice 
corresponds approximately to the inflection point of A^(7r) as 
a function of ti.) In the white area, N{q) has its maximum at 
''"/2 < g < TT, at 7r/2 for large t2. In the exact diagonalization, 
we take Ne — 8 fermions on L = 16 sites. 



less fermions until A^(7r) drops sharply at a level-crossing 
transition at t2 ^ O.lbVf'^^ see also Fig. |31 At the level 
crossing, the ground state changes fundamentally; N{q) 
develops a broad continuum with a maximum between tt 
and tt/2 (moving to tt/2 at large t2) rather than at tt. 
Just as t2 does not influence the charge order, weaken- 
ing with ti, the level-crossing transition driven by t2 is 
hardly affected by ti . This can be seen by comparing the 
ti = 0.02V and ti = 0.071^ curves for spinless fermions in 
Fig. [3l Consequently, the WL phase is bounded by ver- 
tical and horizontal lines in the ti~t2 plane, see Fig.[31 

While the transition between the two CDW phases 
with q = TT and q ^ n depends on both ti and t2, it 
is remarkable that the WL is never affected by the com- 
bination of hopping processes. We would actually ex- 
pect some cooperative effects between ti and t2 because 
NNN hopping is no longer frustrated in the presence of 
ti, see Fig. [T] Due to the DW delocahzation, (b<-^c), 
two-DW states should gain energy with and nonzero 
t2 should thus help destabilize the charge ordering. The 
solution is found in the process shown in (f<->g): For spin- 
less fermions (all arrows in Fig. [1] pointing up) , process 
(a<-S'e) and process (f^g) are equivalent. Since two elec- 
trons swap places in the second case, the resulting Fermi 
sign leads to destructive interference and the lowest-order 
processes associated both with ti and with t2 cancel out. 

After this discussion of the spinless model, we now 
turn to electrons with spin. Due to the dominance of 
the Coulomb repulsion and the classical nature of WL 
ordering, we might not expect charge ordering to be af- 
fected significantly by the spin degree of freedom as long 
as U V . However, the behavior of iV(7r) obtained us- 
ing the DMRG for electrons with spin indicates that there 
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is a surprisingly strong influence even for U = lOOV^. In 
contrast to spinless fermions, where t2 does not affect the 
behavior of A^(7r) as a function of ti, we find the charge 
order to be considerably weakened at ^2 > for electrons 
with spin, see Fig. [2l We can understand this by consid- 
ering the processes of Fig. [T] The states depicted in (c) 
and (g) differ by their sequence of up and down spins. 
Process (b<->c) is then no longer canceled by (f^g) , as it 
is for spinless fermions. Consequently, a kinetic energy 
contribution oc t\t2/f^ is no longer forbidden by the Pauli 
principle. 

Our interpretation is corroborated by analytic consid- 
erations: The additional DW motion due to t2 favors two- 
DW states and changes the gap relevant to Eq. dH from 
A = eo ~ V/2 to A = eo - 2t2. This leads to the dash- 
dotted line in Fig. [21 which indeed describes the weakened 
charge order seen in the DMRG at i2 > 0. For i2 < 0, 
however, the DMRG results are described by the spinless 
gap A = ep. Since spinless fermions are equivalent to the 
fully polarized FM state, this indicates ferromagnetism, 
see below. For U = W and small ti < 0.071^, where AF 
superexchange ~ ^t^/U destroys the polarized state, pro- 
cesses oc t\t2 retain their impact and strengthen charge 
order, see the full line in Fig. [2l 

For fermions with spin, the sharp transition as a func- 
tion of t2 shown in Fig. [3] becomes asymmetric with re- 
spect to the sign of t2- Even for very small ti = 0.02V^, 
the cooperation between ti and t2 is enough to render 
the breakdown of the WL charge order more gradual for 
t2 > than for t2 < 0. For ti = 0.07^, charge order is 
strongly reduced for t2 > 0, and the sharp drop in A^(7r) 
as a function of t2 has disappeared, in stark contrast to 
the spinless model. 



IV. MAGNETISM 

One expects AF interactions for the WL4i due to su- 
perexchange in the generalized Hubbard model. Yet our 
observation that the charge structure factor at negative 
t2 agrees closely with results for spinless fermions, see 
Figs. [2] and [3l is already an indication that the ground 
state in this regime is FM. This is indeed the case, and 
our DMRG studies in fact yield fully polarized ground 
states for some parameter sets at negative t2. As can be 
seen in Fig. [5l the FM interval increases with ti. 

In the following, we analyze the magnetic exchange by 
using perturbation theory valid when ti,t2 V <^ U . 
The robust WL charge order leads to a modulated Heisen- 
berg chain with spins at every second lattice site. Mag- 
netism is therefore described by an effective Heisenberg- 
like Hamiltonian 



7?j = J^(S, -S 



i+2 



1 ^ 

-nini+2) 



(5) 



where i runs only over the even sites, where the L/2 spins 
forming the WL are located. The total ground-state en- 
ergy then is i? = {H,j) -\- EpMi where EpM is the energy 
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FIG. 5: (Color online) Ground state energy Eq for L = 2A 
versus t2 for U = 4V with ti = 0.051^ (circles) and ti = 
0.07V (triangles). Filled symbols indicate a fully polarize 
ground state and horizontal lines the energy of the FM state 
at t2 = 0. The dashed and dash-dotted lines are analytic 
results obtained from perturbation theory, Eq. see text. 



of the fully spin-polarized state, which is equivalent to 
the ground-state energy of spinless fermions. There are 
two distinct mechanisms that contribute to the exchange 
constant J — Jse + Jke- The first term is the usual 
superexchange, which involves a doubly occupied inter- 
mediate state and therefore has the energy scale U in the 
denominator: 



JsE 



4i| 
U 



Utj 8tlt2 



(6) 



The second term Jke — denoted as the kinetic 
exchange — arises from a spin exchange without any 
doubly occupied sites, i.e., exactly from the same effect 
that weakens the charge order for t2 > 0: Quantum 
interference between processes (a<-^e) and (f^g) in Fig.[T] 
is destructive in the polarized FM state and constructive 
in the AF singlet, which leads to an exchange energy 



0^2 

7 ^ 1 
■Jke — — 



1 



1 ~ 2i2/eo 



- 1 



Atlt2 



(7) 



that depends on the sign of t2- The NN correlation 
function of the ID quantum antiferromagnet is given by 
K = -ln2 -I- 1/4 « -0.443, corresponding to (S^ • 81+2) 
for our modulated chain. Inserting k as well as Eqs. © 
and ([7]) into Eq. (O, we obtain an analytic estimate, valid 
at small hopping, for the energy of the AF state 

Eaf ^ ^{JsE + JKE)iK - 1/4) + Efm ■ (8) 

We compare this analytic result to numerical DMRG 
data for ti = 0.05^ and ti = 0.07^, in Fig. [5l This 
figure shows that the analytic curves given by Eq — 
imii{E AF, Epm) closely model the numerical ground- 
state energies for not-too-large t2 with no fitting param- 
eters. Moreover, the analytical boundaries of the FM 
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following from Eq. (|4]) for t2 > have the form 



<J ^-(eo-2t2)v/l/^W-l 



(9) 



FIG. 6: (Color online) Magnetic phases and charge structure 
in the WL regime. FM (o): m > 0.99M; AF (black x): 
m < O.OlAf; 'almost AF' (gray x): O.OIM < m < O.IM, 
where m = 5101(5101 + 1) was obtained by DMRG for L — 
24, U = 4V, and Af = Smax(S'max + 1) = 42. Gray lines give 
potential lines for N{Tr), i.e., the strength of the alternating 
charge order. The lines were interpolated from data points 
obtained at the same parameter values as the magnetic data. 



phase, namely if -it\/U and t\ ^ ~{U/eQ)tl, match 
the magnetic phase boundaries found using the DMRG. 

This may be compared with the boundaries in the 
ti-t2 plane describing the appearance of four Fermi 
points, namely t2 < — (ii/4) sec(7rn/4)^ and ^2 > 
(ti/4) csc(7rn/4)^, where n is the filling. For the quarter- 
filled case, i.e., n = 1/2, only the first relation is of inter- 
est and yields the condition —00 < t2 < — 0.293ti for the 
appearance of FM in the Hubbard model, as observed 
in Refs. and [2^. It is perhaps not surprising that 
the boundary relevant for the Hubbard model, which is 
linear in ti, is completely different from the boundaries 
obtained for the Wigner lattice, which are both quadratic 
in ti. Nevertheless, ferromagnetism at quarter- filling is 
found at negative t2 in both cases! 

The phase diagram in Fig. [6] shows the total spin in 
the ground state of a chain with L = 24 for the range 
of ti,t2 corresponding to the WL regime. The dashed 
lines represent the analytical boundaries t2 and t2 of the 
FM region. Kinetic exchange ([7]) — the only magnetic in- 
teraction surviving for U — > 00 — is AF for <2 > 0. For 
t2 < 0, Jke raises the energy of the AF state over that 
of the FM state, see Fig \5\ In the FM state itself, this 
effective FM exchange is, ironically, absent, because the 
Pauli principle forbids the ring exchange in the polar- 
ized state. At large negative t2 and not-too-large U, AF 
superexchange ^ once again dominates. 

The phase diagram Fig. [S] also contains contour 
lines N{q = TT,ti,t2) — C with the constant C = 
0.9,0.8, ...,0.4, which indicate the strength of the alter- 
nating {q — tt) charge order of the Wigner lattice. In 
accordance with Figs. and [31 we find that the charge 
order dies off most quickly in the singlet states in the 
^2 > region. The analytic contour lines versus ^2 



They agree with the DMRG data, just as Eq. ^ agrees 
with iV(7r) in Fig. [2l Since our analytic result describes 
the unbiased DMRG simulations so well, we conclude 
that the kinetic exchange indeed drives the suppression 
of charge order in the AF regime for t2 > 0. 



V. DISCUSSION AND CONCLUSIONS 

We have investigated charge order and magnetism of 
the ID quarter-filled Wigner lattice with nearest and 
next-nearest neighbor hopping. Starting from the regime 
^ij 1^2! with extremely strong alternating charge or- 
der stabilized by the long-range Coulomb repulsion, we 
find that increasing NN hopping ti drives a crossover to 
a, Akp charge density state with weaker charge order but 
unchanged modulation period, whereas increasing NNN 
hopping t2 leads to a sudden level-crossing transition, de- 
stroying the alternating charge order 4^ For ti, |t2| ^ V 
and spinless fermions, we find that there are no mixed 
processes involving both ti and t2 because destructive 
interference removes the lowest order processes oc iit2- 
Consequently, the WL is bounded by a vertical crossover 
line and horizontal phase transition lines in the ti-t2 
phase diagram for spinless fermions, see Fig. 21 

However, in the case of real electrons with spin, we find 
that processes (x tft2 are absent in the FM state (as for 
spinless fermions) but not in the AF state (as illustrated 
in Fig.[T]). This results in an effective magnetic exchange 
— 4ift2/eQ which favors the AF relative to the FM state 
for positive t2 and disfavors the AF state for negative 
t2- This peculiar effect is corroborated by DMRG data, 
where we indeed find AF as well as FM ground states 
which depend on ti and t2. The phase boundaries of 
the FM phase obtained here for the quarter-filled WL 
are distinct from those obtained for the ti-<2 Hubbard 
models driven by large U and Fermi surface topology, 
where ferromagnetism is found whenever the fully polar- 
ized Fermi sea is split in twoi^ This is perhaps not that 
surprising, as the charge order in the WL is not caused 
by a quantum mechanical Fermi-surface instability but 
by the strong Coulomb repulsion, which does not depend 
on the Fermi surface. Thus, the situation in the WL 
is very different from charge and magnetic order driven 
by Fermi surface instabilities in the Hubbard model with 
purely on-site Coulomb repulsion, which was studied ex- 
tensively in Ref. 16. 

Magnetism in the WL is not driven by Fermi-surface 
instabilities, but actually has more in common with 
the FM reported for a model of coupled chains with a 
symmetry-breaking on-site potential!^ In these coupled 
chains, charge order is stabilized by a strong on-site po- 
tential, which corresponds to the spontaneous symmetry 
breaking by long-range Coulomb repulsion in the case of 
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the WL, and is likewise independent of the Fermi sur- 
face. Starting from strong charge order, we have been 
able to derive the effective magnetic exchange terms us- 
ing perturbation theory, and find the resulting magnetic 
energy to be in good agreement with DMRG data, see 
Fig. [51 The most important terms in the effective ex- 
change Hamiltonian are the AF superexchange oc tj/C/ 
involving a doubly occupied site and the kinetic exchange 
term oc tft2/V^. A related FM exchange mechanism 
based on strong charge ordering and NN hopping has 
been invoked for two-dimensional kagome lattices. 

Since neither charge order nor magnetism are driven 
by a Fermi surface instability but are determined by 
the Coulomb interactions, one might not expect that the 
charge order depends significantly on the magnetic cor- 
relations. Indeed, Wigner lattices are often discussed in 
terms of spinless fermions, and the magnetic exchange is 
added only in terms of a modulated Heisenberg model 
for the given charge order. Our calculations show that 
this picture is too simple: Figures [2] and [3] reveal that 
magnetic correlations have a very strong impact indeed 
on charge order. In fact, the weakening of charge or- 
der in the AF phase at positive t2 is due to the same 
processes oc if ^2 that cause the kinetic exchange. While 
these virtual processes cancel in the FM phase, they in- 
volve domain-wall excitations in the AF phase, see the 
sketch in Fig. [TJ Controlled by the ^2 dependent two-DW 
gap A, the admixture of virtual DW excitations leads to 
a reduction of charge order in the AF state at positive t2 , 
as can be seen in the phase diagram Fig. [SI Within the 
FM regime at negative t2, charge order does not depend 
on ^2- At large negative ^2, however, the antiferromag- 



netic phase reappears, but now, surprisingly, the charge 
order becomes stiffer as \t2\ is increased. This peculiar 
behavior finds its explanation in the t2 dependence of the 
domain-wall gap A in the AF state. 

In summary, we have shown that, in Wigner lat- 
tices with next- nearest- neighbor hopping t2 , a kinetic ex- 
change mechanism is at work, which favors ferromag- 
netism for negative NNN hopping t2 and might ex- 
plain the spin polarization recently observed in strongly 
charge-ordered carbon nanotubeSfi^ In contrast to the 
usual superexchange processes, which involve virtual ex- 
citations across the Mott-Hubbard gap ^ U , kinetic ex- 
change arises from virtual transitions across the Wigner 
lattice charge gap on V U . Although Fermi sur- 
face topology describes neither charge nor magnetic or- 
der in the strongly correlated WL, we nevertheless find 
that quantum interference of electrons is important in 
the WL. In fact, the spin degrees of freedom have such 
a strong impact in the AF regime that charge ordering 
cannot be described reliably in terms of spinless fermions, 
even in the extreme WL regime ii, |i2| ^ V ^ U . As 
this effect is intimately related to the kinetic exchange 
mechanism^ it may also be relevant in higher dimensions, 
e.g., the above-mentioned kagome systems. 
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